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Abstract 

We present chronocoulometry experiments and equilibrium Monte Carlo simula- 
tions for the electrosorption of Br and CI on Ag(100) single-crystal electrode surfaces. 
Two different methods are used to calculate the long-range part of the adsorbate- 
adsorbate interactions. The first method is a truncated-sum approach, while the second 
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deeply. 
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is a mean-field-enhanced truncated-sum approach. To compare the two methods, the 
resulting isotherms are fit to experimental adsorption isotherms, assuming both a 
constant electrosorption valency 7 and also a coverage-dependent 7. While a constant 
7 fits the Br/Ag(100) well, a coverage-dependent or potential-dependent 7 is needed 
for Cl/Ag(100). 

Keywords: Bromine electrosorption; Chlorine electrosorption; Continuous phase transition; 
Chronocoulometry; Lattice-gas model; Monte Carlo simulation. 



1 Introduction 

Specifically adsorbed anions strongly influence the structure and dynamics of adsorbed layers 
on electrode surfaces |2J. Among these systems, halides electrosorbed on single- crystal 
metal surfaces have been extensively studied both experimentally and theoretically over 
the last decade. Examples of these studies include the electrosorption of CI, Br, I, and 
F on low-index surfaces of Ag || |, BJ employing classical electrochemical techniques 
El B, [7J, (see also literature cited in ref. and structure sensitive techniques such as 
electroreflectance |S|], in-situ scanning tunneling microscopy (STM) |], [TTJ, surface X-ray 
scattering (SXS) j|, EJ, [TTJ, and X-ray absorption fine structure (XAFS) 



Ocko, Wang, and Wandlowski || showed that the sharp peak observed in cyclic 
voltamograms of Br electrosorption on Ag(100) corresponds to a continuous phase transition 
in the layer of adsorbed Br. This transition separates a low-coverage disordered phase at 
more negative electrode potentials from a c(2 x 2) ordered phase at more positive potentials. 
Moreover, in recent static and dynamic studies, Mitchell, Brown, and Rikvold [TBI |T3 



have numerically investigated the phase ordering and disordering mechanisms in cyclic- 
voltammetry (CV) and sudden potential-step experiments. Using Monte Carlo simulations 
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of a lattice-gas model first used by Koper |TR |TBJ, they also produced theoretical adsorption 
isotherms for this system that were fit to experimental adsorption isotherms. From these fits, 
parameters such as the electrosorption valency, 7, and the next-nearest neighbor interaction 
parameter, nnn , were estimated. 

Less attention has been given to CI electrosorption on Ag(100), perhaps because previous 
studies || have assumed that this system is similar to Br/Ag(100). In the present paper we 
compare Monte Carlo simulation results with experimental equilibrium isotherms obtained 
by chronocoulometry for both CI and Br, using a lattice-gas model and two different methods 
for calculating the long-range part of the configuration energies. 

The remainder of this paper is organized as follows. The experimental procedures and 
results are discussed in Sec. 0. The lattice-gas model is presented in Sec. [| Equilibrium 
simulations using the two different methods of calculation are discussed in Sec. |], and the 
fitting procedures and resulting parameter estimates for each system are discussed in Sec. ^. 
Finally, a summary and conclusions are presented in Sec. [| 

2 Experimental 

The electrochemical measurements were performed with Ag(100) single-crystal electrodes 
(4 mm diameter and 4 mm thickness) using the so-called hanging meniscus technique. The 
silver electrodes were chemically etched in cyanide solution before each experiment until 
the surface appeared shiny [I7j and, after careful rinsing in Milli-Q water, annealed in a 



hydrogen flame for about 30 s. After cooling in a stream of argon the electrode was quickly 
transferred into the electrochemical cell. Contact with the electrolyte was established under 
strict potential control, usually at values close to the potential of zero charge of 0.05 M 
KCIO4 (E pzc « —0.900 V) The counter electrode was a platinum wire, and a saturated 
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calomel electrode (SCE) in a side compartment served as reference. All potentials are quoted 
with respect to the SCE. The temperature was (20 ± 1)°C. 

The set-up and procedures for the electrochemical experiments (cyclic voltammetry, 
capacitance measurements, chronocoulometry) were described previously We focus 

here only on some additional details of the chronocoulometric experiments : The potential 
was initially held at a value E\ between —1.375 V and —0.300 V (25 mV interval length), 
up to 300 s for the lowest chloride concentrations, and then stepped to the final potential 
Ef = —1.400 V, where chloride is completely desorbed from the electrode surface. The 
waiting time at E\ was chosen (in control experiments) to be always sufficient to establish 
adsorption equilibrium. A comprehensive account of the setup and detailed results of the 



SXS experiments for Ag(100)/Cl will be given elsewhere [ZOfl 



Results of the electrochemical experiments 

Figures 0(A) and|l|(B) show voltammograms (CV) and capacitance (CE) curves (5 mV peak- 
to-peak amplitude, 18 Hz ac-frequency) for the Ag(100) electrode in 0.05 M KCIO4 and 
0.02 M KC1 + 0.03 M KCIO4 as obtained after continuous cycling of the potential between 
—1.40 V and —0.10 V. The capacitance curves merge at E < —1.20 V, indicating complete 
desorption of chloride. The onset of chloride adsorption is characterized by broad peaks in 
the CV and CE curves, which are la belled PI in Fig. [I]. In- situ SXS experiments revealed 
that chloride adsorbs in this potential region as a lattice gas corresponding to randomly 
adsorbed species on fourfold hollow sites P0| . The sharp peak P2 in Figs. [TJ(A) and 0(B) 



represents the continuous transformation of the lattice gas into an ordered c(2 x 2) chloride 



adlayer [ 20| . This phase is stable until rather complex AgCl x -surface compounds are being 



formed at more positive potentials [E > —0.20 V for the system shown in Fig. WA)] |7j, ^] 
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Equilibrium data on the adsorption of chloride ions at Ag(100) electrodes were obtained 
from chronocoulometric measurements in a gently stirred solution. Figure |I](A) shows 
selected charge density curves for KC1 in (0.05 — x) M KCIO4 + x M KC1. Adsorption 
of chloride causes a positive charge to flow to the metal side of the interface. This is 
accompanied by a negative shift of the potential of zero charge. Phenomenologically, the 
charge density vs. potential curves as measured in the presence of chloride concentrations 
Cq > 10~ 3 M are composed of two segments. The region of low charge densities corresponds 
to the potential range of the broad features in the CV and CE curves. The second one 
appears at a potential just positive of the adlayer phase transition at P2, corresponding 
to charge densities between 35 and 40 /zC cm -2 . It is characterized by an initial increase 
of 5q M /5E. Subsequently the slope decreases and seems to become less dependent on the 
chloride concentration for Cq > 10~ 3 M at the most positive potentials studied. 

The corresponding Gibbs surface excess has been calculated at constant potential (see 
Fig. 0) and constant charge (not shown, c.f. P0|), using the approach of Stolberg and 



Lipkowski |2(| 2T|| . The surface excess increases in the potential region of disordered 
chloride and at constant concentration monotonously with potential. At Cq > 10~ 3 M 
one observes in the potential region of the sharp volumetric peak P2 a distinct increase in 
slope (r ~ 6 x 10 -10 mol cm" 2 ). At the most positive potentials experimentally accessible, 
the surface excess levels off and approaches a chloride-concentration independent limit of 
r m = 10.2 x 10~ 10 mol cm -2 . This chronocoulometric result is close to 9.94 x 10~ 10 mol cm -2 
expected for a complete c (2 x 2) chloride adlayer, and is well supported by in-situ SXS 
experiments employing the so-called interference scattering technique |f2T)f . 

The electrosorption valency, as determined from a plot of the charge density vs. surface 
excess at constant potential for the Ag(100)/Cl system is shown as an inset in Fig. [| Using 
this method, the electrosorption valency is found to be potential dependent, decreases with 
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increasing potential, and approaches -0.45 at -0.40 V. This value is larger than the result for 
Ag(100)/Br (7 ~ —0.80) || [L3j (see also discussion in Sec. ||), indicating the stronger ionic 
character of adsorbed chloride, in comparison to bromide on silver. 



3 Lattice-gas Model 

The model used for both the Br and CI systems is a lattice-gas model similar to that used 
by Koper JL5[ [RJ and Mitchell et al. []i~3 , 14|. In this model the Br or CI ions adsorb at the 



four-fold hollow sites of the Ag(100) surface, as shown in Fig. [3]. The model is defined by 
the grand-canonical effective Hamiltonian 

L 2 

H = -^2<PijCiCj ~~pJ2 c i (!) 

i<j i=l 

where Yn<j * s a sum over an pairs of sites, 4>ij are the lateral interaction energies between 
particles on the ith and jth sites measured in meV/pair, and fx is the electrochemical 
potential measured in meV/atom. Here the local occupation variables Cj can take the values 
or 1, depending on whether site i is occupied by an ion (1) or solvated (0). 

The simulations were performed on an L x L square lattice, using periodic boundary 
conditions to reduce finite-size effects. The interaction constants 0y between ions on sites % 
and j a distance apart (measured in Ag(100) lattice spacing units, a = 2.889 A H) are 
given by 



-OO Tij 



''J 



r 

where the infinite value for = 1 indicates nearest-neighbor exclusion, the interaction 
for large is most likely a combination of dipole-dipole and surface-mediated elastic 
interactions, and negative values of <fiij denote repulsion. The previous studies by Koper jl5 



H)| have shown that large but finite nearest-neighbor repulsion has only minor effects on the 



coverage isotherms. Consequently, we have chosen a simplified model with nearest-neighbor 
exclusion. 

In the weak-solution approximation, the electrochemical potential ~p is related to the bulk 
ionic concentration C and the electrode potential E (measured in mV) as 

C 

Jl = -p + k B T In — - e-yE (3) 

where /Z is an arbitrary constant, Co is a reference concentration (here taken to be 1 mM), 
e is the elementary charge unit, 7 is the electrosorption valency |22 |, and /Z has the sign 
convention that ~p > favors adsorption. 
The anion coverage is defined as 

N 

e = N- l Y.c t , (4) 

i=l 

where N = L 2 is the total number of lattice sites. The coverage can be experimentally 
obtained by standard electrochemical methods as well as from the integer-order peaks in 
surface X-ray scattering (SXS) data [|, ^0 . 



To compute the isotherms 6(jt) or 0(E), the change in energy must be calculated before 
each attempted Monte Carlo move. This requires evaluation of the interaction sum in Eq. ([I]), 
which is computationally intensive even for modest system sizes. In the next section we 
compare a direct summation approach with a more efficient truncated summation combined 
with a mean-field approximation for the influence of adparticle pairs at large separations. 



4 Monte Carlo Methods 

The Monte Carlo simulation proceeds as follows. On a square lattice with L = 64 and 
periodic boundary conditions we randomly select a lattice site, i, and attempt to change 
the occupation variable (q = 1 — > = or q = — * q = 1). The energy difference, 
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AH, between the intial state of the system and the proposed state is then calculated using 
Eq. ([[]). We use the Metropolis acceptance probability |f23[ 



TZ = min 



1, exp 



AH \ 



(5) 



for accepting the new trial configuration as the next configuration in the Monte Carlo 
sequence. 

Since the value of AH. determines TZ in Eq. (|5[), the approximations made to calculate 
the large-r^ contributions to the pair sum in Eq. (|I|) are important. In our computations, 
we used two different methods to calculate ATi. The first was a direct-summation-plus- 
truncation method previously used by Mitchell et al. |T3| , |T4|| . In the second method, which 



we introduce here, we include the exact contribution of fewer sites on the surface, using a 
mean-field approximation to calculate the contribution from all the sites at larger distances. 

4.1 Direct summation 

The first method involves calculating the contribution of the sites that are less than or equal 
to five lattice spacings away from the site of interest. Since nearest-neighbor adsorption is 
prohibited, the acceptance probability vanishes for any move that would bring two particles 
into nearest-neighbor positions. In addition to the four nearest-neighbor sites, there are a 
total of 76 sites such that 1 < < 5. These are divided between 12 different distance 
classes as shown in Fig. f|. This summation provides approximately 87% of the total energy 



due to the long-range interactions [14 ]. 



This method of approximation proved to be acceptable for simulating Br adsorption on 
Ag(100) [13|, [14| as it produced values for the electrosorption valency 7 that are in good 



agreement with experimental data. On the other hand, because of its computationally 



S 



intensive nature, this method of calculating A7i made it harder for dynamic simulations to 
achieve CV scan rates within the experimental range |14 . 



4.2 Mean-field-enhanced method 

The method of calculating ATi. can be modified from the previous one so that, without 
reducing the accuracy of the calculation, the cutoff radius is reduced from five to three 
lattice constants. This decreases the number of sites beyond nearest neighbors that are 
included in the summation from 76 to 24, as shown in Fig. ^. 

This modification includes a mean-field approximation to estimate the interaction energy 
contribution from those particles that are farther than three lattice constants from the site 
in question. Consequently, the change in energy associated with desorbing from the surface 
a particle at lattice site i is 



AHi = 2 3 / 2 0r 



Krij<Z'ij \ l<ry<3 ' ij 



(6) 



which is both more efficient and more accurate than calculating the exact energy contribu- 
tions out to Tij = 5. Here Hoc = Y^^Lx 1/r 3 is a Zucker sum over the sites of a square lattice 



of unit lattice constant |15| , and the quantity in the parenthesis is a constant, independent of 
the adsorbate configuration. It can easily be calculated to arbitrary numerical precision. The 
energy change upon adsorption at a previously empty site i with four empty nearest-neighbor 
sites is the negative of the result in Eq. (||). 

This modification of the calculation of the energy changes for the Monte Carlo steps 
makes the program run more than twice as fast as when using the method described in 



Sec. |4J] and used in Refs. |L3|, [14| . It is consequently better suited for dynamical simulations 
fl2l. 
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The results of the calculations performed with this method are similar to those using 
exact summation for r < 5. At higher coverages, the numerical effect of introducing the 
mean-field approximation is more pronounced, as expected. We note that this method is 
successful because the interactions in our model are repulsive, which ensures that the spatial 
distribution of adparticles is approximately uniform on large length scales. For systems with 
attractive interactions more sophisticated methods to calculate the contributions from the 
long-range interactions, such as the Fast Fourier Transform or Fast Multipole Method . 
would be needed. 

4.3 Comparison 

A comparison of the calculated coverage as a function of the electrochemical potential, 
obtained with both mean-field and non-mean-field methods for different cutoff radii, is 
shown in Fig. [5| The mean-field enhancement makes a significant contribution to the 
adsorption isotherms for both cutoff lengths. However, when the mean-field enhancement 
is used, increasing the cutoff radius beyond three does not significantly improve the results. 
Tests with even smaller cutoff radii revealed that three is an optimal choice to increase the 
computational speed with only a minimal loss of accuracy. 



5 Parameter Estimation 

To estimate the parameters in the lattice-gas model, nnn and 7, we performed standard 



equilibrium MC simulations [23] with AH calculated by the mean-field-enhanced method 
at room temperature (k B T = 25 meV or T rs 290 K) to obtain 9(]£) for different 
parameter values. These simulated isotherms were then compared with the experimental 
chronocoulometry data using two fitting procedures. The first procedure assumed a constant 
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electrosorption valency 7, while the second assumed a coverage-dependent 7(0). The 
experimental coverage data for Br/Ag(100) H, which were also used in Refs. flT3| , P|, were 
provided to the authors of those papers by J.X. Wang. 



5.1 Constant 7 

To extract the values of 7 and nnn that would best fit the experimental data, we used a 
least-squares fitting procedure that minimizes x 2 with respect to three fitting parameters for 
CI and Br electrosorbed on Ag(100) @, 

(fc) 

X 2 (0nnn,7,/Io)= E E ^{E^ ft, 7, </W; C7fc) - 0-*(J#; C fc )] 2 (7) 
fc=l Z=l 



where # sim and # exp are the simulated and experimental coverage values respectively, k is the 
concentration index, N conc is the number of different concentrations used in a fit, and I is 
the point index of the experimental data points. 



This fitting procedure, which provides good fits for Br/Ag(100) in Refs. [|T3j, [TJJ (see 
also Table 1) does not give satisfactory fits for Cl/Ag(100). The assumption that the 
electrosorption valency 7 is independent of the coverage is consistent with the Br results, 
but appears to be invalid for the CI system. 



5.2 6 dependent 7 

Since a coverage independent 7 does not appear to be a valid assumption for Cl/Ag(100), 
we assume that 7 depends linearly on the coverage as 

l(d) = lo + lid (8) 

This adds to the model a fourth parameter, 71, that must be determined by comparison to 
the experimental results. The ^-dependent 7 is allowed to vary in the interval [—1,0], which 
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limits the values that 70 and 71 can assume according to Eq. (H). The range of variation for 

70 is [—1, 0], while it is [—2, 0] for 71 since the latter is multiplied by 8, which can assume a 
maximum value of 0.5. As determined by x 2 [Eq. (0)], the inclusion of a coverage- dependent 
electrosorption valency improves the fits to the adsorption isotherms, shown in Fig. || It also 
gives reasonable values for 7 and 4> nrm . The resulting parameter values, based on fitting CI 
and Br simulations to the experimental isotherms, are shown in Tables 1 and 2, respectively. 

The tables include values for the fitting parameters, 7 C for the constant 7 fits and 70 and 

71 for the coverage- dependent fits. Moreover they include values for nnn and x 2 P er degree 
of freedom, x 2 , for all the fits. Fitting to the experimental data was performed both for the 
truncated-sum Monte Carlo isotherms and for the mean-field-enhanced isotherms. For each 
of these simulated isotherms we performed three kinds of fits. First, the numerical isotherms 
were fit to a single set of experimental data with a certain halide ion concentration (columns 
3-8 in Tables 1 and 2, labeled 1-fit). Second, each numerical isotherm was simultaneously 
fit to two experimental isotherms with ionic concentrations of 10 mM and 20 mM for CI or 
1 mM and 10 mM for Br (columns 9 and 10, labeled 2-fit). Third, numerical isotherms were 
fit to three experimental isotherms simultaneously with the concentrations 6, 10 and 20 mM 
for CI and 0.1, 1.0, and 10.0 mM for Br (columns 11 and 12, labeled 3-fit). 

In fitting two or more isotherms simultaneously, several conclusions may be drawn from 
the tables shown. With one exception, the values of 71 in Table 1 for Br/Ag(100) are zero 
(i.e., |7i| < 0.01) or nearly zero for all of the fits. This supports the conclusion that a 
constant 7 is a good approximation for that system. Moreover, the quality of the fits is 
enhanced by using two or more concentrations simultaneously because this utilizes a greater 
number of experimental data points in one fit and constrains 7 more directly through the 
ln(C/Co) term of Eq. (|3]). For constant 7 the 2-fit and 3-fit values for </> nnn and 7 for 
the mean-field-enhanced method are consistent with those obtained by the non- mean-field 
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method. However, smaller values of x 2 are found for the mean-field-enhanced method than 
for the non-mean-field method, indicating that the mean-field-enhanced isotherms fit slightly 
better to the experimental data than the non-mean-field ones. Comparing the values found 
using the 2-fit and 3- fit methods for mean-field-enhanced isotherms, suggests nnn = —21 ±2 
meV and 7 C = —0.71 ±0.01 as reasonable estimates. These results are consistent with those 
obtained in Refs. |5], [13], |I4 |. 

The results for the fits to the Cl/Ag(100) experiments are given in Table 2, whose layout 
is the same as for Table 1. Notice that 71 assumes non-zero values for most of the fits, and 
that - except for the fit to the 6 mM isotherm - the values of x 2 obtained with variable 7 are 
less than those obtained with the constant-7 fits. This supports the conclusion that, for the 
CI system, a coverage-dependent 7 should be adopted. The exception can be attributed to 
the fact that the 6 mM experimental data points do not cover the whole 9 range, especially in 
the important region near the phase transition. Consequently, we hereafter ignore the 6 mM 
data in estimating parameters for the Cl/Ag(100) system. This leaves the 2-fit parameter 
values for a ^-dependent 7. The mean-field and non-mean-field parameters are similar and, 
as determined by x 2 , provide equally valid descriptions of the experiment. From this we 
conclude that, for Cl/Ag(100), nnn = -14±2 meV and 7 = (-O.66±O.Ol) + (-O.68±O.Ol)0. 
While these results for 70 and 71 are more negative than those obtained from the charge- 
density curves in Sec. the trends of a stronger dependence on coverage or electrode 
potential and a higher (less negative) overall value than for Br/Ag(100) are the same. The 
reasons for the quantitative difference between the estimates based on the two methods to 
determine the electrosorption valency are left for future study. 
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6 Conclusions 



In this paper we have presented a new method for implementing the lattice-gas model for 
the electrosorption of halide ions on a Ag(100) surface by introducing a mean-field-enhanced 
method to calculate the contributions to the configuration energy of widely separated 
adparticle pairs. We fit adsorption isotherms obtained by Monte Carlo simulation to the 
experimental adsorption isotherms based on chronocoulometric data after comparing them 
to the results obtained using the non-mean field method. 

The isotherms from both methods were fit to experimental data for both CI and Br. 
It was shown that it was safe to assume a coverage-independent or potential-dependent 
electrosorption valency for Br on Ag(100) while a coverage-dependent electrosorption valency 
had to be implemented for the CI on Ag(100). 
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Table 1: Fitting parameters for Br adsorption on a Ag(100) single-crystal surface, using isotherms computed with a cutoff 
distance of three lattice constants, and with and without mean-field interactions included. The quantity x 2 is the x 2 per degree 
of freedom of the fit. Here, "1-fit" means individual fits for each of the three different bulk concentrations given; "2-fit" means 
a single, simultaneous fit for the two bulk concentrations given; and "3-fit" means a single, simultaneous fit for the three bulk 
concentrations given. 





1-fit; 0.1 mM, 1 mM, 10 mM 


2-fit; 1 mM, 10 mM 


3-fit; O.lmM, 1 mM 10, mM 


mean-field 


no mean-field 


mean-field 


no mean-field 


mean-field 


no mean-field 


7 = Const. 


7c 


-0.51 


-0.49 


-0.46 


-0.47 


-0.46 


-0.45 


-0.70 


-0.67 


-0.73 


-0.70 


</>nnn (mcV) 


-9 


-10 


-10 


-8 


-10 


-11 


-22 


-23 


-22 


-23 


X 2 X 10 5 


0.564803 


2.14327 


7.34073 


0.680078 


2.56305 


8.09803 


6.89425 


9.46991 


7.31550 


9.76821 


7 = 7o + 7i0 


7o 


-0.54 


-0.49 


-0.46 


-0.53 


-0.46 


-0.45 


-0.70 


-0.66 


-0.72 


-0.68 


7i 


-0.01 


0.00 


0.00 


0.02 


0.00 


0.00 


0.00 


0.04 


-0.01 


-0.19 


^nnn(mcV) 


-10 


-10 


-10 


-11 


-10 


-11 


-22 


-21 


-21 


-16 


X 2 x 10 5 


0.54580 


2.20281 


7.54464 


0.69541 


2.63425 


8.32297 


6.98496 


9.53128 


7.36551 


8.68241 



Table 2: Fitting parameters for CI adsorption on a Ag(100) single-crystal surface, using isotherms computed with a cutoff 
distance of three lattice constants, and with and without mean-field interactions included. The quantities shown and the 
organization of the table are the same as in Table 1. 





1-fit; 6 mM, 10 mM, 20 mM 


2-fit; 10 mM, 20 mM 


3-fit; 6 mM, 10 mM, 20 mM 


mean- field 


no mean-field 


mean-field 


no mean-field 


mean-field 


no mean-field 


7 = Const. 


7c 


-0.29 


-0.35 


-0.40 


-0.28 


-0.35 


-0.46 


-0.57 


-0.54 


-0.44 


-0.44 


</>nnn (mcV) 


-4 


-6 


-8 


-4 


-7 


-14 


-19 


-20 


-12 


-14 


X 2 X 10 5 


4.06844 


5.05880 


11.2260 


4.11751 


4.63378 


10.3544 


11.8302 


11.7546 


13.2892 


13.1223 


7 = 7o + 7i0 


7o 


-0.29 


-0.46 


-0.66 


-0.28 


-0.47 


-0.65 


-0.67 


-0.65 


-0.43 


-0.42 


7i 


-0.00 


-0.36 


-0.68 


-0.00 


-0.36 


-0.67 


-0.66 


-0.70 


0.02 


0.04 


</>nnn (mcV) 


-4 


-7 


-12 


-4 


-9 


-14 


-14 


-14 


-11 


-12 


X 2 x 10 5 


4.17015 


3.35255 


5.61327 


4.22045 


2.74380 


4.28522 


5.93654 


5.41895 


13.2942 


13.1225 
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Figure 1: (A) Cyclic voltamogram (50 mV s" 1 ) for Ag(100)/(0.05 - x) M KC10 4 + x M KC1, 

x=0 (dotted curve) and x=0.02 M (solid curve). (B) Capacitance vs. potential curves 

for the system of Fig. |l|(A), 10 mV s _1 . (C) Charge density vs. potential curves for 

Ag(100)/(0.05 - x) M KC10 4 + x M KC1 (ac in M), as indicated in the figure. 
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Figure 2: Plots for surface excess for Ag(100)/(0.05 - x) M KC10 4 + a; M KC1 (x in M), 
employing the electrode potential as the independent electrical variable for selected chloride 
concentrations. The inset shows the potential-dependent electrosorption valency. 
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Figure 3: Br or CI (bigger spheres) adsorbed at the 4-fold hollow sites of the (100) surface 
of Ag (smaller spheres). The squares of the grid frame correspond to the adsorption sites. 
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Figure 4: Sites included in the energy calculation. The dashed lines surround the 76 sites 
included in the direct summation for 1 < r < 5. The heavy solid lines enclose the 24 sites 
included in the direct sumation for 1 < r < 3, which is used with the mean-field-enhanced 
method. 
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Figure 5: Simulated coverage isotherms vs. electrochemical potential, computed with 
different methods. Here R is the cutoff radius up to which exact summation was performed. 
The simulations were performed at T xs 27° C with nnn = —26 meV. The line segments 
connecting the data points are included as a guide to the eye. 
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Figure 6: One simulated Cl/Ag(100) isotherm fit to two experimental isotherms with 
different bulk concentrations. The simulated isotherms were obtained using the mean- 
field-enhanced method described in Sec. 4.2. (a) Using coverage independent 7. (b) Using 
coverage-dependent 7(0) = 70 + 7i#- 24 



